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Four-body interaction plays an important role in many-body systems, and it can exhibit interest- 
ing phase transition behaviors. Historically it was the need to efficiently simulate quantum systems 
that lead the idea of a quantum computer. In this Letter, we report the experimental demonstration 
of a four-body interaction in a four- qubit nuclear magnetic resonance quantum information proces- 
sor. The strongly modulating pulse is used to implement spin selective excitation. The results show 
a good agreement between theory and experiment. 
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I. INTRODUCTION 



Quantum computers have advantages over the classical 
counterparts in simulating quantum systems [l[ and solv- 
ing some hard problems, such as factoring large number 
and searching unsorted databases 01- Among the vari- 
ous candidates for implementing the large- scale quantum 
computer in the future and demonstrating quantum al- 
gorithms to corroborate existing theories, liquid nuclear 
magnetic resonance (NMR) has been proven a convenient 
and practical method to learn lessons for the other physi- 
cal systemsp . NMR quantum computer is an Ising-type 
computer j4| where two-body interactions take the form 
of TTJijalal/2, known as J— couplings in NMR, where cr* 
denotes the Pauli matrix of the i-th spin, and Jij denotes 
the strength of coupling between two spins. 

There have been much interests in many-body inter- 
actions. Besides two-body interactions, many-body in- 
teractions are valuable sources for quantum information 
processors. For example, the three-spin interactions can 
speed up the quantum state transfer in the Heisenberg 
spin chain Four-body interactions have attracted 
much interests recently _6]. The systems with many- 
body interactions can exhibit interesting phase transi- 
tion behaviors Q, such as quantum entanglement phase 
transitions The three spin-interactions in the spin 
chain can induce the quantum criticality that cannot be 
measured by concurrence, because three-spin interactions 
generate three-qubit entanglement Four-body inter- 
action may play an important role in phase transition in 
some condensed matters [6]. 

Simulation of quantum system is one of main applica- 
tions of future quantum computers. Practical factoring 
and searching applications of quantum computer usually 
require hundreds even thousands of qubits. But the simu- 
lation of quantum systems may require only a few dozens 



of qubits. Thus simulating quantum systems may well 
be the first practical application of the early practical 
quantum computer. It is helpful now to study the sim- 
ulations of quantum system with a few qubit quantum 
information processor, to locate problems and gather ex- 
periences, in particular the unitary operations and the 
extent of decoherence in existing apparatus. In fact, it 
has been successfully demonstrated that three-body in- 
teraction can be simulated very well in NMR quantum 
computers 0, It is a reasonable assumption that the 
four-spin interactions relate to four-qubit entanglement, 
which is still unclear for us currently. In this paper we fo- 
cus on implementing the four-spin interactions in a NMR 
quantum computer. Our work is a valuable step in ex- 
ploring quantum simulations, and also a crucial step for 
implementing the quantum phase transitions induced by 
four-body interactions and investigating the relation be- 
tween the four-spin interactions and four-qubit entangle- 
ment experimentally. 



II. GENERATING FOUR-SPIN INTERACTIONS 
USING NMR 



Our task is to decompose the four- spin evolution into a 
series of one-spin operations and J-couplings. The one- 
spin operations are realized by radio frequency pulses. 
The J-coupling evolution 



-i{Tr/2)JijTal(Tl 



(1) 



can be realized by standard NMR spin-echo techniques 
i. Under U,,{T), at evolutes as 
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where 9 — TrJijt. 



Through some calculations [10| one finds that the n- 
spin interaction can be decomposed as the {n — 1)- spin 
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interactions by iteration 
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By introducing 
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Eq. ([3|) can be further expressed as 
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From Eq. Q one finds that the many- spin interaction 
can be decomposed into the operations that can be di- 
rectly reahzed by NMR. When n = 4 one obtains 

= g~i(7r/4)<7^g-i(7r/4)a>^g-i(7r/4)ajg-i(,r/4)a;j 
X e-'('r/4)T>Jg-i(7r/4)<T;;g-i(7r/2)Ji234Ta>^ 
-i(7r/4)<T;; -i(7r/4)o-^<7^ i(7r/2)o-;; i(7r/4)<7;^ 
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where J1234 is the effective strength of the four- spin in- 
teraction. The above equation can also be represented 
as 



X e-*('^/4)T>^g-»(7r/4)aJg-i(,r/2)Ji234T<T>^gi(7r/4)<Tj 
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III. IMPLEMENTATION 



We use Carbon-13 labelled crotonic acid dissolved in 
D2O as the sample. The chemical sketch of crotonic acid 
is shown as Fig. [U where CI - C4 are assigned as qubits 
1-4, respectively. The protons are decoupled during the 
whole experiment. The experiments are implemented on 



a Bruker DRX 500 MHz spectrometer. The temperature 
is controlled at 22 °C. 




FIG. 1; (Color online) Molecule of crotonic acid. The blue 
spheres denote protons and the green spheres denote oxygens. 
The chemical shifts are ui =21468.9 Hz, 1/2 = 15255.6 Hz, 
V3 = 18668.0 Hz, and V4 = 2190.4 Hz. The J coupling con- 
stants are J12 =72.4 Hz, J13 =-1.3 Hz, J14 =7.0 Hz, J23 =70.3 
Hz, J24 =-1.6 Hz, and J34 =41.3 Hz. 



The Hamiltonian of the NMR system reads 
4 



(9) 



fc=i 



k<l 



where - are the resonance frequencies of CI - C4. 
The coupled-spin evolution between two spins is denoted 
as 



\Tkl\ 



(10) 



where k/l — 1, 2, 3, 4, and k ^ I. [tm] can be realized by 
averaging the coupling constants other than Jki to zero 
0. The pulse sequence to implement [1/2 J12] is shown 
in Figl21 where the evolution time 1/2 J12 is divided into 
eight identical segments. 
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FIG. 2: Refocusing scheme to implement [I/2J12]. The el- 
lipses represent the strongly modulating tt pulses. The evolu- 
tion time 1/2 J12 is divided into eight identical segments. 



We use Eq. ([7]) to implement Uzzzz{T) through the 
pulse sequence 
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23 
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(11) 



'2Ji2 



where [f denotes a tt/2 pulse along — y axis on C2. 

The corresponding evolution is e^^^^^^^'^y . All spin selec- 
tive pulses are strongly modulated pulses (SMPs) 
A SMP consists of a series of non-selective (hard) pulses 
that modulate the system's dynamics strongly to produce 
precisely a desired spin-selective unitary propagator. In 
our experiments the fidelity of each SMP is larger than 
0.99. The total duration time of the whole experiment is 
about 80 ms. 

We choose the state 

P^n^ = ^1, (12) 

as the initial state, which is prepared by 

^ [gradl - [|]^ 

from the thermal equilibrium [l^ . Here we use the 
deviation density matrix to describe the state of the NMR 
system [l^, [l3| • The carbon spectrum for the system in 
Pini is shown in Fig. [3l where the signals are chosen as 
the reference signals for the following spectra. The small 
J13 causes the partial overlapping peaks. 



18740 18720 18700 18680 18660 18640 18620 18600 Hz 
FIG. 3; The C3 spectra when the system lies in pini — at. 



Under the four-body interaction, the state of this sys- 
tem changes from pini to 



al cos(7r Ji234r) -f alala^a^^ sa\{TT J i23iT) . 



Pzzzz (^) 

(13) 

When 7rJi234r — mr/A with n = 0, 1, 8, the spectra of 
C3 are shown as Fig. |4l The experimental results agree 
on the theoretical expectations. By integrating over the 
eight peaks we obtain the evolution of {a"^) as a func- 
tion of time, 7rJi234r, and is shown in Fig. [S] The 
curve can be fitted as (cr^) = 1.081 cos(1.0087rJi234T), 
which agrees well with the theoretical expectation (cr^) = 
cos(7r J1234T). The small discrepancy is due to the imper- 
fection of pulse and decoherence. 
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FIG. 4: NMR spectra of C3 in evolution under a four-body 
interaction. From top to bottom, the figures correspond to 
7rJi234T = nn/A with n = 0, 1, 8, respectively. 



IV. SUMMARY 

We have experimentally simulated the four-body inter- 
action in a four-qubit NMR quantum information proces- 
sor. The experiment results show good agreement with 
the theoretical expectations. The SMP makes the sim- 
ulation in NMR very well. With this experiment, one 
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FIG. 5: The expectation value (a^) versus ■KJ1234T 
under XJzzzz- The curve can be fitted as (a^) = 
1.081 cos(1.0087rJi234T). 



can proceed to demonstrate large quantum system sim- 
ulations, and look into interesting physical phenomenon 
such as phase transitions in quantum systems with four- 
body interaction. The simulation method used here in 
NMR techniques can be generalized to other Ising type 
quantum computer. 
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